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ABSTRACT 

The kinematics of stars and planetary nebulae in early type galaxies provide vital clues 
to the enigmatic physics of their dark matter halos. We fit published data for fourteen 
such galaxies using a spherical, self-gravitating model with two components: (1) a 
Sersic stellar profile fixed according to photometric parameters, and (2) a polytropic 
dark matter halo that conforms consistently to the shared gravitational potential. The 
polytropic equation of state can describe extended theories of dark matter involving 
self-interaction, non-extensive thermostatistics, or boson condensation (in a classical 
limit). In such models, the fiat-cored mass profiles widely observed in disc galaxies are 
due to innate dark physics, regardless of any baryonic agitation. One of the natural 
parameters of this scenario is the number of effective thermal degrees of freedom 
of dark matter {F^) which is proportional to the dark heat capacity. By default we 
assume a cosmic ratio of baryonic and dark mass. Non-Scrsic kinematic ideosyncrasies 
and possible non-sphericity thwart fitting in some cases. In all fourteen galaxies the 
fit with a polytropic dark halo improves or at least gives similar fits to the velocity 
dispersion profile, compared to a stars-only model. The good halo fits usually prefer 
values from six to eight. This range complements the recently inferred limit of 
7 <^d < 10 (Saxton & Wu), derived from constraints on galaxy cluster core radii 
and black hole masses. However a degeneracy remains: radial orbital anisotropy or 
a depleted dark mass fraction could shift our models' preference towards lower F^; 
whereas a loss of baryons would favour higher F^ . 

Key words: dark matter — galaxies: elliptical — galaxies: haloes — galaxies: 
individual (NGC 821, NGC 1023, NGC 1344, NGC 1400, NGC 1407, NGC 3377, 
NGC 3379, NGC 3608, NGC 4374, NGC 4494, NGC 4564, NGC 4697, NGC 5128, 
NGC 5846) — galaxies: structure — stellar dynamics 



1 INTRODUCTION 



Simulations that ; 
matter successfully 



ssume simple, 
reprodu ce 



non-interacting dark 
^ cosmic large-scale struc- 
ture (see e.g. ISpringel et all 120061 '). However, they also 
predict power-law cusps of dark density at the cen- 
tres of all virialised structures (e.g. Dubinski fc Carlbere 
1991'; 'N avarro et "al]ll996l : iMoore et al.ll 19981 : iDiemand et al. 
200^, .2OO5I '). Dynamical and other circumstantial evidence 



Blumcnthal ct al. "l986l;'Gncdin 
2004; Roinano-Di'az ct al 



for 



of 



uniform 



dark density 


(e.g Flores & Primackl 


19941: Moord 


1994 


Burkerd Il995 


: IWeldrake et all l2003l: 


Gentile et al.l 


2004 


Simon ct al.' 


2005: de Blok 20051: iKuzio de Narav et all 


2006: Wciimans et al.ll2008l). It remains debatable whether 



Zhao"2002':'Gncdin ct alj 
2008; Jardcl & Scllwood 2009). 
The observational situation in elliptical galaxies is less 
clear than for other morphological types, because of the rel- 
ative difficulty of observing suitable kinematic tracers at 
a range of radii. These data may suffer degeneracies be- 
tween mass p rofiles and the orbital a nisotropics of kine- 
matic tracers (|Binnev fc MamonI Il982l ). X-ray emission (if 
detected) is informative if hydrostasis is a fair assump- 
tion. Observations of some X-ray emitting ellip tical 



the present-day cores 
olent, primordial form 
how conspired to 



are innate, or whether a vi- 
of baryonic "feedback" some- 
erase (or exacerbate) all cusps (e.g. 
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favou r the presence of massiv e dark halos (e.g. Awaki et al.l 
19941 ; iLoewenstein fc White! 1 19991 ; iMathews fc Brighentil 
2OO3I ; iHumphrev et al.ll2006h . For some controversial galax- 
ies, dark matter appears dynamically absent, perhaps more 
diffuse than in standard cosmogonies, or possibly disguised 
by stellar orbital anisotropics (e.g . Mendez et al.l 2001 ; 
Rornanowskv et al.ll20ol : lDekel et al.ll2005l ; lde Lorenzi et al l 
20091). In other cases the in f erence of dark matter is cleare r 
(e.g lo'Sullivan fc PonmanI |2004| : Ide Lorenzi et al.l |2008| V 
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Gravitational lens models suggest steep, isothermal slopes 
around the observed radii. However, inner regions show 
moderate dark mat t er densities (e.g. Ger hard ct alj I2OOII: 

2006 : Thomas et all 



Ferreras et alJ 20051: Cappellari et alJ . 

2OO7I : Weiimans et al. 2009t ). despite the cuspy pred 



tions of simple coUisionless dark matter cosmogonies, and 
the expectation of adiabatic contraction due to the in- 
fall of the baryons to wards the centre of the halo (e.g. 
iBlumenthal et al.lll986l '). 



Before N-body simulations became popular and highly 
resolved, flat-cored proflles like the "pseudo-isothermal 
sphere" were by default assumed in studies of galaxies 
and galaxy clusters. Observations of galaxy-scale dark cores 
now motivate the examination of alternative dark matter 
theories where the physics naturally cause cores, regard- 
less of any baryon i c "feedback" (e.g. Spergel fc Steinhardtl 
2000l:|Peeblesll2000l: lArbev et al.ll20q3l:lAhn fc Shapirdl2005l: 



Saxton fc Wul l2008i r TNunez et al.l (|2006l ') and IZavala et al] 
( 20061 ) compared the scaling relations of disc galaxies using 
eithe r a standard dark matter profile (NFW, iNavarro et al.l 
1 19961 ). or a polytrope - intended to represen t coUisionless 
dark matter subject to Tsallis thermostatistics (|Tsallij|l988l ) 
- and found the l atter to explain better the observations. 
iBohmer fc Harkd (|2007l ) also found good fits to the rotation 
curves of low surface brightness galaxies and dwarf galaxies 
with a polytropic dark matter halo (which they justify via 
the formation of a Bose-Einstein condensate). 



Furthermore, polytropic dark matter halos can equally 
represent other interesting scenarios, such as self-interacting 
dark matter, where the core has formed due to the sup- 
port of dark pressure in the dark field domain. Polytropic 
models are naturally consistent, self-gravitating distribu- 
tions of classical mass, and may be seen as analogues to 
adiabatic fluids. Some early simulations with numerical ap- 
proximations to d ark hvdrodvnami cs (iMoore ct al.' 20_0; 
lYoshida et "al]l200Q : .Burk crt _2000,; Kochanek fc White 20~ 
entailed (dark) thermal conduction leading to gravothermal 
catastrophe, and hence isothermal cusps that were sharper 
than those of coUisionless simulations. However it was even- 
tually shown (iBalberg fc Shapird '200^: 'B alberg et alll2002l : 
lAhn fc Shapirdl2005l ) that the simulations were unwittingly 
(and inappropriately) initiated in near-collapse conditions. 
It was shown that conductivity is only significant for SIDM 
of medium interactivity. Weaker SIDM has infrequent colli- 
sions and conduction is slow. Strongly self-interactive dark 
matter has vanishing conductivity (as it returns to an adia- 
batic limit) which defers collapse well beyond the age of the 



In this paper we e xtend the polytropic dark halo model 
of lSaxton fc Wul (|2008l ) to galaxy scales. We apply the model 
to fitting planetary nebula and stellar kinematics of ellipti- 
cal galaxies. We present fits for a number of cases, including 
"stars only" models, "stars + halo" models with isotropic 
orbits and a constant global baryon fraction, and finally 
we also consider the issue of orbital anisotropy and a non- 
universal baryon fraction. 



2 METHOD 

2.1 Spherical model and its solutions 

We assume a fixed stellar mass distribution, wit h a projected 
surface density approximately described by the lSersig (|l968l ) 
profile. Sersic-like profiles appear ubiquitous among ob- 
served galaxy spheroids, although the reasons why this form 
emerges in nature have not yet been demonstrated analyti- 
cally. For the three - dimen sional stellar density, we adopt the 
Prugniel fc SimienI lll997[) spherical deprojection (see also 
Lima Neto et al.ll 19991 : ICiotti fc Bertinlll999l : lMarquez et al.1 
20001 ). 



Pir(x) = paX ^exp ^—b (^x^^" — 1 



(1) 



with dimensionless coordinate x = r/ra, and {b,p) are con- 
stants depending on the index n. The scale radius (3D) is 
approximately equal to the observational, projected (2D) 
half-light radius, ~ Re- The stellar mass within some ra- 
dius is given by 



m*(< x) = 47rn6"fP-^'eVcr;? r[n(3 - p),bx^^"] 



(2) 



involving lower incomplete gamma functions. 

Within this empirical stellar distribution, we introduce 
a dark matter halo, with a density distribution to be solved 
consistently with the global gravitational potential. The lo- 
cal equation of state is polytropic, 



(3) 



with a pressure Pd, density pd, pseudo-entropy Sd and adi- 
abatic index 7d = 1 -|- 2/Fd. This is the natural behaviour 
expected for an isotropic medium comprising classical par- 
ticles with _Fd effective degrees of freedom. Such a con- 
dition is expected if dark matter is self-interacting (e.g . 



Spergel fc Steinhardt 


I2OOOI: Peeblesll2000l: Arbev et al.ll2003l: 


Ahn fc Shapird 12005 


: ISaxton fc Wull2008l). An equivalent 



polytropic equation emerges if the dark hal o is governed 



riasri 



by Tsallis' non-extensive thermostatistics (e.g . iTsallisI 19881: 
Plastino fc Plastind 1 19931 : IZavala et all I2OO6I : iNunez et al] 



If —2 < _Fd ^ 10 then the dark halo may truncate at 
finite mass and radius (Rh)- For any particular model, we 
find this radius by assuming dark hydrostasis 

dPd _ Gm Pd 
dr ' 
where we define the interior mass m — m^{r) -\- md{r) for 
brevity. Computationally, it is more convenient to deal with 
the dark velocity dispersion (or temperature) o"! = Pd/ Pd, 
which obeys 



(4) 



dr 



Fd + 2 



Gm 



(5) 



We take Sd as global constant, and numerically integrate 
the dark profile from the origin outwards. We assume non- 
singular conditions at the origin, 



dpd 
dr 







(6) 



precluding the complication of a compact central mass and 
its sphere of infiuence. At some point the dark velocity dis- 
persion becomes small (also pd — >■ if Fd > 0) and the inte- 
grator switches to use ad as the independent variable instead 
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of r. Integration proceeds exactly to the point where = 0, 
i.e. the natural truncation radius _Rh. We call this the "dark 
surface". The low- mass region between the dark core (say 
where the log-slope of density is —2) and the dark surface 
shall be called the "halo fringe." 

Once the total dark mass Md and Tih are known, we can 
obtain the external gravitational potential, stellar pressure, 
potential energy and thermal energy profiles by integrating 
a Jeans equation (and associated ODEs for the mass profile 
and potential) inwards from infinity. If we define a stellar 
pressure in terms of the radial velocity dispersion of stars, 
P^^ = p*o-J then we have the spherical Jeans equation 

dPj, _ Gmp^ 2/3P» 

dr r 

where /3 = 1 — erg /crj is the local stellar velocity anisotropy. 
When integrating ((Tjl far outside the dark halo, we change 
the independent variable to u = 1/r. Having obtained all 
the relevant constants of integration at the dark surface i?h, 
we can easily evaluate profiles of the stellar and dark com- 
ponents in both the interior and exterior. 

2.2 Model specification and projection 

In any theory with dark matter degrees of freedom 
fixed, the pseudo-entropy Sd and central dark tempera- 
ture are formally free parameters. When discussing or fit- 
ting families of physically related models, it is more useful 
to prescribe the stellar mass fraction within some radius, 
= ^* I (j^-k + 'Tid ) , or the dark fraction /id = 1 — /i* . 
If galaxies retain all their primordial endowments of matter 
and darkness, then total stellar value would approximate the 
cosmic baryon fraction /i^f(oo) ~ 0.1(2 When seeking par- 
ticular numerical solutions, we shall attain a chosen [i^ioo) 
by fixing the dark entropy Sd and adjusting the central (j\ 
iteratively. Then (for a given stellar background) there is an 
implicit relation between Sd and the dark truncation radius 

Model observables are calculated by integrating suitable 
3D stellar properties at radii r = + where h is the 
projected radius on the sky plane and z is the line-of-sight 
coordinate. The column density of stars is simply 



= 2 



and the sightline integrated pressure is 

n* = 2 / (1-/3) 



dz 



(8) 



(9) 



The mean line-of-sight velocity dispersion is then ctios ~ 
\/II^^/E^. The orbital anisotropy is crucial to the stellar 
observables. We will mainly consider models with perfect 
isotropy (/3 = everywhere). In later sections ( i)3.5p . we 
will consider anisotropy varying radially according to the 
Osipkov-Merritt model (|Osipkovll 19791 : lMerrittlll985h . 



/3(r) = 



f.i _|_ di 



(10) 



^ We note that the latest WMAP-5 r esults give a best fit for a 
cosmological baryon fraction of ~0. 2 jPunklev et al.ll2009l'l . We 
refer the reader to §§3.4 for an estimate of a change in this num- 
ber. 



with the scale radius a treated as a fitting parameter. 

In practice, integration of the stellar profile to radii 
at large orders of magnitude eventually samples a tenu- 
ous, hot region where drops much faster than P^. For- 
mally, the integrated velocity dispersion rises greatly. Phys- 
ically, this contribution should not occur, because the den- 
sities imply fewer than one star within the relevant volume. 
Therefore we cut off the integrals (jSj and ((9} at some three- 
dimensional radius far outside the kinematic tracers of a 
particular galaxy. 

Calculations are performed in natural units of the Sersic 
profile: pe = Re ~ G = 1. A numerical minimal-x^ fitting 
procedure automatically yields the normalisation for a par- 
ticular data set. This auto-normalisation counts as an extra 
free parameter, decrementing the "degrees of freedom" in 
the reduced scores. 



3 RESULTS 

We fit the kinematics of stars and planet ary nebulae of 
a sam ple comprising twelve galaxies from ICoccato et al.l 
(|2009t ) as well as NGC 1400 and NGC 1407. These data 
are separated into major and minor axis components. We 
fit the axes separately, to obtain an indication of the appli- 
cability of our model's spherical assumption. Additionally, 
we fit kinematic data for two of the galaxies without de- 
composition along axes: fo r NGC 821 we extract data from 



iRomanowskv et al. 



from lDouglas et al 



120031 ). For NGC 3379 we use figure 11 
I2007h. The data for NGC 1400 and 1407 



are taken from Proctor et al.l |200^. Notice in this case the 
data are extracted over a circular aperture, so only one axis 
can be considered. 

We assume that the stel l ar ma ss component follows the 
profile of lPrugniel fc SimienI (|l997h . with the Sersic param- 
eters shown in the left three columns of Table [T] For nine of 
the galaxies, these para meters come from ph otometric fits 
to Sersic-only models bv lCoccato et all (|2009|. their table 3) . 
For NGC 3379 we take n = 4.7 from lDouglas et all (|2007li . 
For NGC 5128 we assu me n = 4 and obtain _Re from the B- 
band photometric fit of lDufour eral] (Il979l). For NGC 5846 , 
we convert Sersic parameters from iMahdavi et al.l (|2005h . 
For all g alaxies, we keep the d istances and luminosities tabu- 
lated bv lCoccato et alT(|2009l ). Initially we attempt fits with- 
out a polytrope ("stars only") to test the necessity of dark 
matter, and to provide benchmarks for the fitting qualities 
of halo models. 



3.1 Fits witii stars only 

Previous modeling of PN kinematics have suggested a 
poverty of d ark matter i n the halos of som e isolated ellipti - 
cal galaxies (|Mendez et al. 2001: Ro manowskv et al.l |2003| ) 
or that the halos are less concentrated than popular coUi- 
sionless dark matter cosmogonies predict (e.g.lPouglas et al.l 
I2OO7I : iNapolitano et al.ll2009l : ICoccato et al.l l2009l ). To pro- 
vide baselines for comparison, we first fit the Coccato et al. 
data with a purely Sersic model (using the published, photo- 
metrically obtained parameters) and no dark matter. Except 
for NGC 4697, the goodness of fit scores - given by the re- 
duced Xr = x'^/'^ ^ w-^ll above one, and therefore bad 
fits (see columns 4 and 5 of Table [1] and the solid lines in 
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figure [T]). The model profiles of the projected velocity dis- 
persion are smoothly declining outside a half-light radius, 
whereas the observational profiles may be flat, or they may 
dip and rise with radius (Figure [T]). Some of these complex- 
ities and departures from ideal Sersic mass proflles may be 
responsible for the formally poor fits. 



3.1.1 Improved stellar fits 

Using the photome tric Sersic parameters tabulated by 
ICoccato et ahl (|2009t ). some galaxies give much worse val- 
ues for the reduced than others. In order to improve 
these cases, we try fitting isotropic Sersic models without 
dark matter, and adapt the parameters (n, _Rc) for a better 
kinematic fit, (n', R'^). In this search, we omit the planetary 
nebulae, relying on absorption line kinematics only. 

Table[T]gives the alternative fits for axes 1 and 2 of each 
galaxy, befor e and after kinematic fitting. The left columns 
are based on ICoccato et"al] (|2009l ') photometric parameters; 
the right columns (primed quantities) are our kinematic fits. 
We give the quality of fit (Xr){i,2} and the total (stellar) 
mass m^k-^j ^ji where the subindex refers to the axis consid- 
ered in the analysis. These fits are shown in figure[T]as dotted 
lines. In some cases we use subindex to refer to data that 
were obtained within circular annuli. We use an "amoeba" 
algorithm to find better fitting parameters, {n',R'^), with 
considerably improved (x?)' on each axis. For NGC 3377, 
1407 and 5128 the amoeba runs away forever along a degen- 
erate x^ valley, heading towards i?c — >■ oo. 

In the majority of attempts at (n', R'^) kinematic fits, we 
find Rc drifting to large values beyond the data, resulting in 
unphysically large total stellar masses. NGC 4564 and 3608 
have smaller R'^ than the photometric values. NGC 1023 is 
the only galaxy where the photometric and kinematic fitted 
Re values are moderate and comparable. Disagreements of 
photometric and kinematic Sersic parameters might indicate 
spatially varying stellar mass/light ratio, or the existence of 
dark matter. Alternatively, the disagreements may be in- 
fiuenced by degeneracies between the parameters, and the 
incomplete, finite radial range of the data. 



3.2 Fitted halos with isotropic orbits 

For the ICoccato et al] (|2009l ) galaxies, data are available 
along two principal axes. We fit these independently, using 
a stellar mass profile with Sersic model parameters adopted 
from the published stellar photometry (i.e. n and Rc from 
table [l]). The central properties are tuned so that the dark 
halo and stars are globally in a chosen ratio. We explore 
cases with fixed Fd, and vary the halo radius (via Sd) to 
achieve the best fit. Comparing the fits for the two axes 
gives some indication of the quality of the spherical approx- 
imation. (It is an imprecise comparison however, since the 
minor axis data tend to span a shorter radial range, and are 
less restrictive.) For six galaxies there are formally satisfac- 
tory fits, with Xr < 1 (i-e. in NGC 1407, 4374, 4564, 4697, 
5128 and 5846). For the rest, the minimal Xr never drops be- 
low a few (formally unfavourable). This may be partly due 
to ideosyncratic stellar kinematic substructures or noisiness. 

To address this possibility, we repeat the kinematic fit- 
ting process with the stellar aios{r) profile smoothed to a 



best-fit cubic function in ln(r). PN data are left untouched. 
Smoothing the stellar kinematics does not shift the min- 
ima significantly, but lowers the attainable Xr- Thus the 
minima are sharper and deeper, but the landscape is un- 
changed at medium and high levels. Results for the 
smoothed fitting are shown in Figures [5] and [3] The colours 
from violet to red represent dark degrees of freedom Fd = 
{1, 2, 3, ...7, 8, 9, 9.5, 9.9}. The PN velocities were assumed to 
be isotropic. We assume a universal cosmic fraction of dark 
and baryonic mass (i.e. st ars for our galaxies) according to 
the WMAP5 cosmology (|Dunklev et al.ll2009t ). In figure [3] 
we show four galaxies for which we have data within circu- 
lar apertures ( Romanowskv et "al]|2003l : iDouglas et al.|[2007l : 
IProctor et aLlboogl ). 

The usual x^ curve of a galaxy has two minima: a shal- 
low dip around _Rh ~ 0.1-Ro to 1-Re, and a deeper dip at 
radii outside the observations. The physical interpretation 
of the inner dip implies a dark halo more compact than 
the stellar distribution. This would be surprising in the con- 
ventional cosmogonies in which radiative cooling made the 
star-forming gas more concentrated than the dark compo- 
nent. Numerically, at inner-dip configurations the "halo" is 
taking over the role of the stars influencing the inner aios 
proflle. The outer minima seem more physically plausible. 
Sometimes the outer dip is effectively a lower limit on Rh 
(e.g. for NGC 4494) but in other cases the fitting prefers 
a finite range (e.g. less than lOOJJo for NGC 821). If the 
data are clear and widespread enough, the radial variation 
of the slope of pd in the observed range is sufficient to distin- 
guish _Fd and extrapolate the truncation radius. For Fd < 9 
the x^(^h) curves tend towards the "stars only" level as 
-Rh — >■ oo. A large- _Rh halo (of fixed fid) is spread thinly and 
has negligible effect on observable kinematics. The conver- 
gence is closer for lower Fd, as these models resemble the 
incompressible limit (_Fd ~ 0) which has uniform density 
within r < _Rh. The highly compressible -fd 9 cases be- 
have quite differently, as discussed in ^3.3\ 

Figure|4]shows the two-dimensional 90, 95 and 99% con- 
fidence levels for six galaxies in our sample, in the parameter 
space spanned by the dark matter degrees of freedom (Fd) 
and the halo truncation radius (Rh). The crosses mark the 
best fit, and the arrow represents the region over which the 
velocity dispersions are fit. The above mentioned degener- 
acy seen in the ID plots of figure [2] is actually a continuous 
one between Fd and i?h. 

For three exceptional galaxies (NGC 1023, 3608, 4564) 
the inner minimum is very significant, or even dominant, in 
terms of ■ Inspection of the ctios proflles shows that the 
photometrically-derived Sersic model flts them especially 
poorly. The inner part of the stellar profile is unusually cen- 
trally peaked. NGC 3377 also has a peculiarly shaped aios 
proflle, which may be responsible for the slight preference 
for Fd — 9.9 on the first axis, and for the failure to improve 
on the stellar-only baseline fit to the second axis. The mod- 
elling of these four misfit galaxies implies very small stellar 
mass to light ratios (in the B-band), typically T:^ < O.lTf 
(see NGC 1023(i,2), NGC 3377(i) and NGC 4564{i,2) in Ta- 
ble [2]) . This is implausible for the aged stellar populations 
of early- type galaxies. We deprecate the discussion of these 
particular galaxies. 

Except for these odd cases, the "stars + halo" fits tend 
to reject Fd > 9. Where the data are clear and smooth 
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enough, the isotropic modelhng generaUy prefers the dark 
halo to have _fd = 8 or fewer effective thermal degrees of 
freedom. There is a mild preference for the top end of the 
range. In models with -Fk 9 the central dark density is 
insensitive to variations of the outer halo radius (-Rh) or the 
halo specific entropy (sd). This inflexibility in fitting central 
conditions may cause the poor fits for Fd ^ 9. The best fits 
favour halo radii of at least tens of Rc. 

We now comment on the favoured isotropic models for 
specific galaxies. 



3.2.1 NGC 821 

The kinematics of planetary nebulae aro und NGC 821 show 
a rad ially declining velocity dispersion (|Romanowskv et al.l 
This has been taken as evidence disfavouring the pres- 
ence of a dark halo (unless the PN orbital anisotropy has 
significant radial variations) or else the halo is less concen- 
trated than CDM expec tations. 

On both axes of the lCoccato et al.l (|2009l ') data, we find 
that polytropic halo models fit significantly better than stars 
alone. The improvement is best on the first axis. The best 
fits (Xr ~ 1-7) are for _Fd ~ 8 and an outer radius 60 < 
Rh/Re ^ 90. These minima improve significantly upon the 
"stars only" fit, in the case of axis 1 by up to Axr ~25. 

The upper left panel of Figure [3] presents the equivalent 
reduce d-x^ scores obtained by fi tting polytropic dark halos 
to the lRomanowskv et all (|2003l ) data for NGC 821. These 
data were annular in distribution (not divi ded onto axes) 
and less detailed than ICoccato et al.l (|2009l l. The fits also 
prefer ~ 8, but the minima are shallower. For _Fd = 8 the 
fit improves gradually for larger _Rh, but for Fd ^ 7 there is 
an optimum at several tens of Rc. For the dark mass fraction 
within IR e, the most lik e ly val ue is 0.2 — 0.4, compara- 
ble to the ICoccato et all (|2009l ') value (si 0.4). The best fit 
stellar mass-to-light ratio T:^ ~ 5Tq is typical of metal- 
rich p opulations of age 5-6 Gyr (see e.g. iBruzual fc CharlotI 
I2OO3I V 



3.2.2 NGC 1023 

Our analysis favours a dark halo with F^ « 9.9 (near max- 
imum) and small radius Rh < lORc- The poor halo fits for 
this galaxy are perhaps unsurprising given the poor fits to 
the "stars only" basic model. The stellar aios profile is more 
centrally pea ked than the (photometric) Sersic model. It 
was observed l|Noordermeer et al.ll200 j : ICoccato et al I2OO9I ) 
that the inner regions r < 100") are affected by strong ro- 
tation. NGC 1023 is a SBO galaxy with a fast rotating bar 
which would imply a m aximal (i.e. baryon dominated) disc 
l|Debattista et al.H200j '). We omit NGC 1023 in the analysis. 



3.2.3 NGC 1344 

Fits to data from axis 2 are bad mainly because of the 
fiat plateau in velocity dispersion in the inner region. 
Axis 1 gives an acceptable fit for F^ ~6, with a signifi- 
cant improvemen t over the stars-only model. According to 
iTeodorescu et al.1 (120051 ). there is a 3.8x lO^M© dark matter 
halo out to S.5Re. A rough estimate of the baryon fraction - 
assuming 90% of the stellar mass is contained within 3.5-Ro 



(|Graham et al.ll2005l ) - gives a baryon fraction within this 
region of 0.29, in agreement with our estimates (/s ~ 0.3). 

3.2.4 NGC 1400 

This lenticular galaxy is one of the tw o bright mem bers of a 
nearby group (along with NGC 1 407, lGouldlll993l). W e ex- 
tract the stellar kinematics from iProctor et al. I (|2009l ). We 
take their assumed distance, but use th e B magnitude and 
Sersic model from ISpolaor et al.l (|2008l ). Assuming a cos- 
mic baryon fraction, none of the polytropic halo models fits 
as well as "stars only." Perhaps this is attributable to the 
discrepancies be tween gipg profiles f rom different observa- 
tions (noted by iProctor et al.l l2009t ) . We also notice that 
NGC 1400 has a significa nt amount of rotation {v/a ~ 0.5 
at Rc , iBertin et al. Il994l ) , which could possibly explain its 
-t-0.3 dex residual with r espect to the Fundamental Plane 
(|Prugniel fc Simien|[l996l ). 



3.2.5 NGC 1407 

This EO galaxy is the brightest member of a nearby group 
(including NG C 1400. [Could 1991) . The stellar data were 
extracted fromjProctor et al.l (fioogt ). and we refer again to 
ISpolaor etal] (|2008l ') for the light profile. The polytropic 
model with cosmic dark matter fraction gives a good fit, 
with _Fd ~ 7 and truncation radius ~15-Re. The halo models 
with Fd > 9 have significantly poorer fits than for low and 
medium effective degrees of freedom. 



3.2.6 NGC 3377 

Stellar data for this fiattened (E5) galaxy show a pronounced 
"S-bend" in the crioa profile. All the fits are formally bad. 
On the first axis, a halo with high Fd brings a slight - but 
perhaps insignificant - improvement. On the second axis, a 
dark halo does not improve the fit at all, and the Fd = 9.5 
and 9.9 curves are worse than the others. We suspect that 
complicated orbital structure (especially the undulating pro- 
jected kinematics) renders straightforward spherical models 
inadequate and inappliable to NGC 3377. Details of this 
profile have been attributed to the influence of a centra l 
massive object and discy dynam i cs (iGebhardt et al.l l2003l : 
ICopin et aLll2004h . lCoccato et all (|2009l ') note a twist in the 
PN kinematics. We omit NGC 3377 from consideration in 
our conclusions. 



3.2.7 NGC 3379 

This galaxy is a n other of the notable systems that 
iRomanowskv et al.l (|2003l) observed to have a radially de- 
clining (Tios profile. For the lCoccato et al] (|2009l ) data on two 
axes, the dark halo does not improve the flts over a Sersic- 
only model. Dark matter is allowed if the halo has a large 
radius and low central concentration. Cases with Fd ^ 9 are 
much worse fits than for lower Fd. 

For comparison with the split-axial results, we also 
fit the projected velocity dispersion data representing 
circular rings , whic h are extracted from figure 11 of 
Douglas et al.l (l2007h (who in corporated stellar data from 



Statler fc Smecker-Hanelll999l ). A purely stellar model fits 
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poorly: if the orbits are isotropic then Xr ~ 4.1. If dark 
matter is present in the cosmic fraction, then it improves 
the fitting somewhat; the best Xr ~2 occurs for Fd ~ 7 (see 
lower middle panel, Figure [2}. Isotropic models with ^ 9 
fit unacceptably worse than those with Fd ^ 8. 

The simple Sersic model of the stellar mass may be 
inadequate to improve upon these fits. A bump in the stellar 
(Tios profile suggests the presence of kinematic substructure. 

3.2.8 NGC 3608 

Including a dark halo does not improve on stellar-only fits. 
The minima - which have very high values of Xr ~ either 
imply a very compact halo {R ~1-Re) with _Fd > 9 or an 
extended halo (~ 500-Re) with Fd = 1. We tentatively blame 
this outcome on the centrally peaky aios profile. If the Sersic 
profile describes the stellar mass distribution poorly, then 
the fitting r outine compensates b y seeking a high central 
dark density. ICoccato et al.l (|2009l ) remark that this galaxy 
had few PN detections, and possible contamination from 
NGC 3607 nearby. We omit NGC 3608 from further detailed 
examination. 



I2OO6I : Ide Lorenzi et all I2OO8I : ICoccato etall 120091 '). Despite 
the non-sphericity, our fits to both axes favour the presence 
of a dark halo with Fd ~ 7 or 8 consistently. The extrapo- 
lated outer radius is large, _Rh ^ 100-Re. Thus the character- 
istics of this galaxy are consistent with the other apparently 
good fits. 



3.2.13 NGC 5128 

For the famous nearby galaxy NGC 5128 (Centaurus A) the 
axis 1 data strongly favour the presence of a dark halo with 
cosmically mean composition. The best fit indicates Fd ~6 
and R ~ 25i?c. This is fewer dark thermal degrees of free- 
dom than in our other plausible fits. The halo radius also 
seems more compact than usual. The prominent dust lane 
implies a recent merger, which might have left persistent 
departures from spherical hydrostasis in the dark sector as 
well. (Persistent dark pulsations or rotation could possibly 
change the apparent Fd obtainable from inapplicably static 
kinematics.) Unfortunately the results from axis 2 are in- 
conclusive due to insufficient data. 



3.2.9 NGC 4374 

This El radio galaxy (Messier 84) is a member of the Virgo 
cluster. Fits to this galaxy favour the presence of a dark 
halo, with Fa « 8 and R ^100 to lOOOiJo. Unusually, for the 
second axis, the _fd = 9 case improves on the basic "stars 
only" model (although _Fd « 8 is still the best fit). 



3.2.10 NGC UH 

This is another of the intriguing cases where 

iRomanowskv et al.l (|2003l ) found a declining velocity 

profile i mplying a dark matt er halo with a very low 



de nsitv jNapolitano et al.ll2009l ). From data on both axes 
of ICoccato et al~ I (|2009l ). the fitting characteristics resemble 
those of NGC 3379: dark matter brings no improvement. A 
halo is allowed if its outer radius far exceeds the observed 



region. 



3.2.11 NGC 4564 

Models fitted to this E6 galaxy strongly prefer the presence 
of dark matter, but with radius of only a few R^,. Perhaps 
this surprising conclusion is due to the centrally peaky stel- 
lar (Tios profile. The mass and orbital distributions may de - 
part significantly from a Sersic sphere. ICoccato et al.l (|2009l ) 
indicate that there is a rotation curve resembling "an SO 
galaxy rather than ... an elliptical galaxy". The unphysi- 
cally low values of M/L and of the baryon fraction confirm 
the fit for this galaxy must invoke non-spherical models. We 
omit NGC 4564 from the analysis. 



3.2.12 NGC 4697 

This flattened elliptical (E6) galaxy has especially numer- 
ous PN measurements. Its velocity dispersion declines in 
the outskirts, an d the degree of rotation appears to de- 
cline outside Tie (|Mendez et al.ll200ll . [20oi : ISambhus et all 



3.2.14 NGC 5846 

NGC 5846 is the central galaxy of a nearby group, and it 
ap pears round (morp hological type EO-El) and non-rotating 
( .Coccato et al.|[2009l ). We should have expected decisive fits. 
The Sersic para meters are deba t able: b y fitting stellar and 
PN photometry, ICoccato etakl (120091') find n = 12 ± 2, 
R^ = 2903 ± 192''. However iMahdavi et al.l (|2005l ') fit the 
light profile with n = 3.95 ± 0.05 and Re = 69.5 ± 5.5". We 
adopt the latter parameter pair (high Sersic indices, nor- 
mally paired with very large half-light radii are usually in- 
dicative of a fitting degeneracy) . Our fits to axis 1 marginally 
favour a dark halo with a small outer radius, J?h ~ 7.8Ra, 
and Fd ~ 1. On axis 2, the Sersic-only fit is similar to a 
stellar -|-dark matter model. 



3.3 Features of the dark matter profile 

In a Lanc-Emden sphere unaffected by other mass compo- 
nents, the quasi-entropy can take any positive real value, 
and there is a one-to-one relation between the entropy and 
the natural outer truncation radius of th e halo. We find that 
this is untrue for a halo perturbed by a lPrugniel fc SimienI 
|l993) stellar component. For Fd > 6 or so, we numerically 
find that there is a minimum Sd at which the halo fringe be- 
comes infinitely diffuse, failing to truncate. For many models 
the cosmic baryon fraction is unattainable beyond a partic- 
ular range of the entropy (even if the halo is finite). 

The poor fitting by Fd > 9 halos has a related cause. 
When the degrees of freedom are numerous, the core is 
smaller and denser compared to the outskirts. This dense 
core tends to outweigh the central stellar density, provid- 
ing a poor fit to the inner data. For Fd > 9 the dark mass 
always outweighs the stars within 17?o (see upper lines in 
figure [5]). Furthermore, when Fd ^ 9 we find that inflat- 
ing the halo (extending the fringe radially) does not reduce 
the central dark density significantly. For lower Fd, the core 
expands somewhat as the truncation radius increases, and 
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this enables fits with plausibly low dark densities within the 
half-light radius. 

Figure [6] illustrates the halo density profile for the best- 
fitting models of axis 1 of NGC 821 for the choices of as 
shown in figure [2] (with the same colour coding). In our fits, 
this galaxy possesses a significant dark mass fraction within 
the half-light radius, and the dark density at l_Ro is com- 
parable to the stellar density there. The Fd 9 curves are 
structurally distinct: the dark density follows that of the 
stars closely for several decades in radius. The truncation 
radius is large, but nevertheless the central dark density is 
high. The more favourable fits, with ^ 8, involve ha- 
los that truncate within 100-Re, and there is a much larger 
inner region where pd > p* (e.g. r «0.7i?G for Fd = 8). 
For all fits, the logarithmic slope of the dark matter halo 
density is around —1 around and within the half-light ra- 
dius. This should not be mistaken for the pd oc cusps 
of coUisionless N-body simulations. The innately flat core 
of the polytropic halo has been gravitationally pinched by 
the stellar mass distribution. The total mass density index 
wavers around —2 out to several Rc, which may be consis- 
tent with_jm£li£ations from some gravitational lensing stud- 
ies (|Rusin et al.ll2003l : iKoopmans et al]|2006l : iFerreras et al.l 
120081 ). If -Fd < 9 and if our models have typical mass concen- 
trations then the polytropic halo of a truly isolated elliptical 
galaxy would eventually be found to decline more sharply 
than oc r~^, somewhere beyond ~ lO-Ro- Satellite galaxies 
and correlated cosmological macrostructure in the vicinity 
may confuse the practical detection of this truncation. 

Figure [7] shows the scatter of the baryon fraction with 
respect to stellar masses for the best-fit models that ade- 
quately improve over the stellar-only baseline cases. Solid 
dots represent those galaxies with a good fit {Xr < 3). The 
figure shows that within IRe most galaxies appear baryon 
dominated except for massive ones, mainly NGC 1407. At 
5-Re the baryon fraction decreases in most cases to ~20%. 
For comparis on, the analysis of strong gravitational lenses 
from SLAGS (|Bolton et al.ll2008l ) is shown as a grey shade 
in the left panel - the SDSS-selected lenses mainly probe 
the inner (_R < R^) regions. There is a lso consistency with 
the estimates of lCappellari et al.l (|2006l ) using the SAURON 
sample - they quote a dark matter fr action of ~ 30% in th e 
inner Rc; or with the recent work of iTortora et al.l |200^, 
with similar values. The grey cross and error bars in fig- 
ure[7](rigAt) correspond to the lensing analysis of the galaxy 
B1104-181 (at a redshift z=0.73) from IFerreras et all (|2005h . 
which was observed within R < 3.7 Re (limited by the po- 
sition of the lensed images of the background source, which 
controls the accuracy of the lensing estimates). 

3.4 Stellar brightness and mass fraction 

Though we have only fourteen galaxies (and not all of them 
with satisfactory fits) it is interesting to consider their ap- 
parent properties as a population. Figure |8] shows the re- 
lation between B-band stellar mass-to- light ratio (T.*.) and 
the total stellar masses (m^^), which are obtained jointly 
from the fitting process. Those galaxies with a good fit 
(Xr < 3) are shown as solid dots. In grey we also show the 
dynam ical M/L values from Ivan der Marel fc van DokkumI 
l|2007l ) for the galaxies in our sample. Those values are typ- 
ically obtained within Rc, and agree rather well with our 



estimates using a polytropic dark matter halo. For compar- 
ison with stellar populations, we also show in the right side 
of figure[H]the expected stellar M/L for a couple of synthetic 
stella r pop ulations f r om th e GB07 models (e.g. lBruzual~Al 
[2003) for a lChabrieij (|2003l ') IMF. Two metallicities are con- 
sidered, as labelled. We find that the more massive galaxies 
tend to have older stellar populations with ages compatible 
with more detailed work based on equiv alent widths (see e.g. 
iThomas et alll2005l : iRogers et aLll2009l ) 

Among our standard calculations, NGC 821 shows per- 
haps the most convincing evidence for a dark halo of cosmic 
mass fraction. In a set of exploratory calculations, we vary 
the assumed ratio of stellar to dark mass for this particular 
galaxy. As illustrated in Figure [21 the stellar mass fraction 
p^k.(co) does have some effect. If the galaxy is richer in dark 
matter than the cosmic mean, then the best model shifts 
towards slightly higher F^, e.g. for p.^ = 0.05 - a value 
motivated by the lates t estimates of the baryon fraction in 
galaxy halos (see e.g. iMoster et al]|2009l : [Guo et al.ll2009l ) 
- the optimum is around Fd « 9 with 7?h ~ 457-Re. If the 
galaxy is poor in dark matter then _Fd and -Rh both shift 
to lower values. As p^^ — >■ 1, the curves collapse trivially 
to the poor fit of the "stars only" model. The best x^ in- 
dicates moderate depletion of dark matter, e.g. when fixing 
p* = 0.30 we find xl = 0.98 and i?h ~ i.94Re for F^ = 3. 
This would imply halo truncation within the span of the ob- 
served planetary nebulae. Models with p.^(cxj) < 0.2 do not 
imply such a peculiarity. 



3.5 Effect of radial orbital anisotropy 

To explore the role of orbital anisotropy, we try re-fitting 
some of the galaxies assuming the Osipkov-Merritt profile 
for /3(r) (see equation (TDJ . In this model, stellar orbits are 
isotropic in the centre, a nd increasingly radial in the fringe 
beyond a scale radius a (|Osipkovlll979l : iMerrittI 119851 ). For 
"stars only" fits the OM radial anisotropy profile provides 
negligible improvement if the scale radius a is large (ap- 
proaching isotropic models) , and for small a it makes the x^ 
scores significantly worse. Therefore radial anisotropy in the 
outskirts seems unlikely to cure the apparent dark deficit 
(or low halo concentration) of cases like NGC 4494. In this 
sense, "stars alone" models may be insufficient without non- 
spherical effects. 

Figure [TD] shows the effect of Osipkov-Merritt radial 
orbital anisotropy profiles upon "stars -I- halo" models of 
NGC 821 (axis 1) . NGC 821 is one of the cases with a clear 
need for a dark halo and this persists when the fringe is ra- 
dially anisotropic. The x^ curves still discriminate between 
_Fd sharply. The curves prefer low F^ if a < 2; for a = 2 
the optimum is Fd = 3; for a = 3 it is -fd = 6. Thus the 
Osipkov-Merritt radial anisotropy model favours small Fd if 
the scale radius is comparable to the half-light radius, while 
isotropic models favour Fd « 8. The effects of the OM model 
are similar for other galaxies, such as NGC 3379. Overall x^ 
levels are lower for larger a, and thus the isotropic limit is 
formally favourable. To convincingly constrain the amount 
and functional form of the anisotropy profile, it would be 
necessary to fit more detailed models with axisymmetry or 
fitting velocity kurtosis data. Such complexity is beyond the 
scope of our present paper. 
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3.6 Exotic models: -2 < Fd < 

It is wortli mentioning tlie mathematical existence of a range 
of more exotic polytropes, with negative degrees of freedom. 
Models with —2 < Fd < and negative pressure (Pd < 0, 
Sd < 0) describe a "generalised Chaplygin gas," which has 
been proposed to unify dark matt e r and dark energy (e. g 
iBento et all |2002| : iBilic et all |2002| : iBertolami et al.l l2004h . 
Polytropes with this composition have infinite radius and 
mass, so the y need arbitrary truncation a t some background 
density (e.g. IBertolami fc Paramo3l2005l ). Thus it is unclear 
how to naturally constrain the dark mass fraction in Chap- 
lygin models of galaxies. 

However, if the pressure is positive, then the ma ss 
and radius of the halo are finite (jViala fc Horedtlll974al lbl: 
lKimurall98ll : lLipscombel2008l '). In this "anti-Chaplygin gas" 
halo, density is minimal at the origin, and increases to in- 
finity at the outer radius Rh where temperature drops to 
zero. (Farther out, the dark density stays zero, as it does 
outside normal polytropes with < -Fd < 10.) These pecu- 
liar bubble-halo models are calculable, though the infinitely 
dense surface causes difficulty for the numerical integrators 
that project ctios. We performed a sparse set of exploratory 
calculations in this regime, and find a continuation of the 
trends in {Rb, Fd) already seen at the low end of the 
Fd > domain. For NGC 821, the optimum rises (poorer 
fitting) and continues smoothly shifting to smaller i?h as 
Fd \ -2. 

The limit of Fd — —2 describes an isobaric condition. In 
this case, hydrostasis requires the gravitational field to van- 
ish everywhere. There are no bound, finite, static, spherical, 
adiabatic solutions. The choice of F = —2 and Pd = could 
describe a collisionless medium, but the dark density and en- 
tropy Sd would then be arbitrarily spatially variable (which 
is a complication beyond the scope of this paper). 



4 CONCLUSIONS 

We have fitted published kinematics of fourteen diverse el- 
liptical galaxies, assuming a Sersic stellar mass profile plus 
a spherical polytropic dark matter halo that is distributed 
self-consistently in the shared gravitational field. This non- 
standard model is justified by our current lack of knowledge 
about the nature of the dark matter particle. A non-zero 
cross section for non-gravitational interactions or internal 
degrees of freedom will justify this interpretation. Further- 
more, studies on the effect of gravitation on an ensemble of 
particles suggest a polytropic equation of state is required 
even if we are dealin g with a non-interacting, elementary 
particle (|Tsallislll98i '). 

Our fits imply that more massive galaxies have higher 
stellar mass-to-light ratios, and higher central fractions of 
dark matter. We also impose a constraint on the degrees 
of freedom for the dark matter particle, Fd > 6, in agree- 
ment with a previo us study focussed on galax;y clusters 
(|Saxton fc Wufl2008l ). The favoured Fd increases to higher 
values if elliptical galaxies are poorer in dark matter (and 
the favoured Fd is smaller if baryons are lost). Dark matter 
fractions well above or below the cosmic fractions may be 
explicable in several ways. Dark matter may be lost via ab- 
lation or tidal stripping, especially in dense environments. 



Baryonic matter might be lost due to blowout from stellar 
or AGN feedback, or from inefficient cooling (at the massive 
end), letting the baryons escape as warm or hot IGM/ICM 
unbound to the galaxy. 

The tentative preference for six to eight dark matter 
thermal degrees of freedom is interesting when compared to 
circumstantial evidence from galaxy clusters. In their analy- 
ses of galaxy clusters com prising polytropic dark matter and 
radiative gas. lSaxton fc Wu (2008j) found that the condition 
for stationary solutions implies a lower limit on the central 
mass. Compatibility with the scale of observed supermassive 
black holes (as well as the dark matter core radii of clus- 
ters) implies a constraint of 7 < Fd < 10. (The upper limit 
is necessary to form any finite self-gravitating body.) Taken 
together, the evidence from elliptical galaxies and galaxy 
clusters therefore implies that dark matter has Fd ~ 7 to 8 
(if it is the same substance in both contexts). 

Half of the available galaxies are well fit by the poly- 
tropic halo model, despite its serious physical simplifica- 
tions. Some galaxies may require tailored, non-spherical 
modelling or special treatment of the kinematic relics of 
their ideosyncratic merger histories. We have assumed that 
the adiabatic dark halo is well mixed and is non-rotating: 
this may be untrue for a galaxy that has merged less than a 
few dynamical times ago. Our present modelling precludes 
the possible effects of a central star cluster or giant black 
hole, either of which could induce a parsec-scale spike of 
dark and stellar matter (within some gravitational radius of 
influence). This could raise the central velocity dispersion. 
Allowing this feature to emerge in future models might im- 
prove fits to the three galaxies with centrally spiky aios{r) 
profiles. 
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Table 1. Reduced values for isotropic fits omitting PN data, with a Sersic (stellar) mass componen t only. There a re too 
few data on axis 2 of NGC 5128, so in this case we fit both axes together. NGC 1400 and 1407 data JProctor et al.|[200a'l 
are not split by axes. Primed values correspond to a kinematic fit to the Sersic parameters. 
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2 In these galaxies, the kinematic fit to the Sersic parameters does not converge. 

Table 2. Characteristics of the formally best fits for our stellar + polytropic dark matter model for each galaxy. Columns 
from left to right are: the galaxy identifier (subscripted with the axis - axis meaning circular annuli); the dark matter 
thermal degrees of freedom; quality of fit for stars + halo {Xr)i quality of fit with stars alone (x^ ^; including PNe data); 
truncation radius of the halo; stellar mass to light ratio in B band; baryon fraction of the interior (r < IRe)', baryon fraction 
of middle regions (r < 5i?c). The uncertainties are quoted at the la level. 
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Figure 1. Projected velocity dispersion profiles of tiie isotropic models ttiat best fit the ICoccato et aD l|2009l 'l data (on two principal 
axes, red and blue). Stellar data (O) determine these "stars only" fits; the PN (x) data are shown for reference. Solid lines are "stars 
only" models based on the standard photometric Scrsic parameters; dotted lines reflect parameters adapted to fit the kinematics better. 
The dashed lines are the best fits of the model including a polytropic dark matter halo. On the horizontal axis, the observed effective 
radius and the one obtained from the kinematic fit are shown as Rc and R'^, respectively (for some galaxies the "kinematic" Rc extends 
beyond the range of the figure, and for three cases (NGC 1407; 3377 and 5128) the kinematic method does not converge; see text for 
details). 
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Figure 2. Quality of isotropic fits to the Coccato et al. kinematic data for two principal axes, as a function of halo outer radius i?^. 
Colours represent cases of different degrees of freedom, = 1,2,3,4,5,6,7,8,9,9.5,9.9 as indicated in the legend in figure [Sj Dotted 
lines indicate the quality of fit for the Sersic only model. Arrows show the radial coverage of kinematic data. Stellar data have been 
softened to a best-fit cubic. Unusually, some of the minima of NGC 1023, 3608 and 4564 imply very compact halos (i?h $ ^e) fitting 
better than the "stars only" level. More plausible fits (lower two rows) prefer a halo radius outside the observed stars and PN. 
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Figure 3. Same as figure [2l for data without axial splitting llRomanowskv et al.ll20o3 : iDouglas et al]|2007l : iProctor et al.ll2009h . 
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log(R,/Re) 

Figure 4. The confidence levels in the parameter space spanned by truncation radius (-Rh; horizontal) and number of degrees of freedom 
of the dark matter equation of state (Fj; vertical) are showrn for six of the galaxies from our sample. They correspond to the 90, 95 and 
99% (thick like) level, using the likelihood defined by the fit to the velocity dispersion data. The horizontal arrow delimits the range 
over which the observed data are available. The crosses in each panel give the best fit. Notice the best fit for the truncation radius of 
NGC 4697 falls outside of the range shown (see table [2J. 
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Figure 5. The dark mass fraction within li^e, as function of halo outer radius, for three galaxies with adequate fits. Open and closed 
circles mark the best fits to axis 1 and 2 respectively. For _Fd = 9, 9.5 and 9.9 the dark core outweighs the local stars no matter how far 
out the halo extends (i£h)- The horizontal dotted line indicates the regions of baryon-dark matter equality. 
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Figure 6. The best-fit 3D halo density profiles for NGC 821 (axis 1) data. The first panel shows the pd{i") curves shaded as usual according 
to their values, The dashed black line traces the stellar density profile for comparison. The density normalisation po = p^(l-Ro) is 
obtained in the optimisation of each fit. (The t otal masses differ be tween cases, but the stellar mass fraction is fixed.) The dotted black 
line corresponds to a standard NFW profile jNavarro et al.|[l996f l. The second panel shows a quantity dmd/rf(lnr), which emphasises 
the layers where most of the dark mass resides. Visually, the total mass is proportional to the area under each curve. The F^ ^ 9 fits 
concentrate much of the dark mass within the half-light radius, while < 9 fits concentrate dark mass farther out than several Re- The 
NFW model implies that most of the halo mass resides at the largest radii. We also include a choice of negative F^ (thin line), which 
results in an increasing outward density (at decreasing temperature). 




Figure 7. The baryon fraction - defined as m^/{m^ + 'm^) is shown within IRe (left) and 5Re versus total stellar mass. We show 
the best fitting cases that improve over the corresponding "stars only" models. Error bars are shown at the 68% level. The sol i d dots 
correspond to galaxies with a good fit Xr < 3. For comparison, the shaded area in the left panel is the result of iBolton et al. I 1I2OO8I ') 
for gravitational lenses in the SLAGS survey. The grey cross with 95% error bars is lens B1104-181 from [p^rreras et al whose 
lensing and stellar masses are obtained within 3.7-Rc. 



Polytropic dark halos of elliptical galaxies 17 



K o _ 

V ~ 



too 



-1 1 — I — r- 



— 1 1 1 1 r 



Ze 2.5Ze 

^-5 



i-iH 



_l I I L 



J 



_I I I 1_ 



o 

--ID 



10" 

Total Stellar Mass (M^) 



Figure 8. The B-band stellar M/L ratio is shown versus total stellar mass. The points correspond to the best fitting cases that improve 
over the corresponding "stars only" models (er ror bars shown at the 68% confidenc e level). The solid dots are those galaxies with a 
Xr < 3. The grey crosses are M/L values from Ivan der Marel fc van DokkumI | |2007| ) for those galaxies that overlap with our sample. 
The arrows on t he rig ht give the stell ar M/L values for a Simple Stellar Population from the CB07 population synthesis models (e.g. 
iBruzual AlboO'lt) for a lChabrieJ l l2003f) IMF. Two metallicities are considered, as labelled. 
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Figure 9. Reduced curves for halo models fitted to axis 1 of NGC 821, but with different stellar mass fractions, fj,^{oo). 
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Figure 10. Reduced scores of model fits to NGC821 axis 1, with curves of different .F^ coloured as in Figure[2] Each panel shows a 
different choice of the Osipkov - Merritt anisotropy scale radius a (in units of Re). 



